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Nucleation and growth mechanisms of short chains of carbon atoms on single-layer, hexagonal 
boron nitride (h-BN), and short BN chains on graphene are investigated using first-principles plane 
wave calculations. Our analysis starts with the adsorption of a single carbon ad- atom and examines 
its migrations. Once a C2 nucleates on h-BN, the insertion of each additional carbon at its close 
proximity causes a short segment of carbon atomic chain to grow by one atom at at a time in a quaint 
way: The existing chain leaves its initial position and subsequently is attached from its bottom end 
to the top of the carbon ad-atom. The electronic, magnetic and structural properties of these chains 
vertically adsorbed to h-BN depend on the number of carbon atoms in the chain, such that they 
exhibit an even-odd disparity. An individual carbon chain can also modify the electronic structure 
with localized states in the wide band gap of h-BN. As a reverse situation we examined the growth 
of short BN atomic chains on graphene, which attribute diverse properties depending on whether B 
or N is the atom bound to the substrate. These results together with ab-initio molecular dynamics 
simulations of the growth process reveal the interesting self-assembly behavior of the grown chains. 
Furthermore, we find that these atomic chains enhance the chemical activity of h-BN and graphene 
sheets by creating active sites for the bonding of various ad-atoms and can act as pillars between two 
and multiple sheets of these honeycomb structures leaving wider spacing between them to achieve 
high capacity storage of specific molecules. 

PACS numbers: 68.65.-k, 73.22.-f, 75.75.-c, 81.07.-b 



I. INTRODUCTION 

After the synthesis of grapheneP^ a monolayer of 
sj9 2 -bonded honeycomb structure of carbon atoms, in- 
terest has focused on two dimensional nanostructures 
having honeycomb structures. Planar hexagonal boron 
nitride (h-BN) is an ionic honeycomb structure consisting 
of alternatively bonded boron and nitrogen atoms and is 
an analog of graphene. 3 4 Despite the structural similar- 
ity, h-BN differs from graphene with its wide band gap 
and dielectric properties!^ Various boron nitride struc- 
tures like nanosheetsp' nanotubes^ and nanowires^ have 
already been synthesized. Also, recent studies show that 
h-BN can be used to improve the current voltage prop- 
erties of graphene transistors by improving the mobility 
of electrons in graphene as compared to graphene films 
on silicon substratesP These properties hold promise for 
novel technological applications of h-BN structures. 

A thorough understanding of the properties of h-BN 
structure and investigating its functionalization is impor- 
tant for future applications of this wide band gap mate- 
rial. Carbon, being in the same row of the periodic table 
with boron and nitrogen, is one of those foreign atoms 
which can greatly change the physical and chemical prop- 
erties of h-BN. In a recent study™ the effects of ad- atoms 
adsorbed on BN were investigated using first-principles 
calculations, and it was shown that high coverage of car- 
bon ad-atoms can change the magnetic properties and 
band gap of the system, whereas ad-atoms induce local- 
ized states in the band gap at low coverage. 

Single carbon ad-atom and mono-atomic chains of car- 



bon atom are interesting entities, which can function- 
alize h-BN. A short segment of carbon atomic chain 
containing n C atoms, indicated as CAC(n) hereafter, 
is a truly one-dimensional carbon allotrope, and has 
recently drawn attention due to its linear geometry, 
high strength, size-dependent quantum ballistic con- 
ductance and interesting electronic properties. These 
properties of CACs were both theoreticalljM^ and 
experimentall}Sm^ investigated. T heoretica l studies re- 
vealed the Peierls distortion in C ACsp^ 13 * 19 * half- metallic 
and spintronic prope rties j 1 6 " 1 8 1 size dependent quan- 
tum conductance^^ and other geometrical structures.^ 
Chain structures of other group IV elements and group 
III-V compounds were also treated.^ Concomitantly, 
carbon atomic chains are synthesized pDHH Recently, it 
was also shown that CACs can be grown on graphene 
and modify its properties,^ which have also been justi- 
fied by the images obtained earlier using high resolution 
TEMpSl 

In this paper we studied the growth of short carbon 
and BN chains on single layer h-BN and graphene, re- 
spectively. We first examined the adsorption of single 
carbon atom on h-BN by calculating its energy landscape 
and diffusion barrier. This is followed by the investi- 
gation of the nucleation and growth processes of CACs 
on h-BN. We performed both conjugate gradient calcu- 
lations and molecular dynamics simulations in order to 
determine the stabilities and bonding properties of these 
CACs and show how they can grow on the plane of BN 
as new carbon atoms are introduced by one atom at a 
time at the close proximity of an existing CAC. Once 
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FIG. 1: Energy variation of single carbon atom adsorbed on 
various sites of single layer, 2D hexagonal boron nitride struc- 
ture (h-BN) calculated in (4 x 4) supercell. (a) h-BN honey- 
comb structure on which the adsorption energies are calcu- 
lated. Nitrogen and boron atoms are represented by blue and 
green balls, respectively. The most favorable binding site of 
carbon ad-atom is marked by the red star in the figure. The 
path of diffusion of carbon ad-atom with the minimum energy 
barrier of ~ 0.65eV is indicated by stars, (b) Complete en- 
ergy landscape of C ad- atom on h-BN structure. Light blue 
regions show favorable sites and the energy barrier further in- 
creases as the color goes to dark blue and purple, (c) Energy 
variation of carbon ad-atom is shown along the path indicated 
by red arrows in (a). The energy difference between the most 
favorable site (indicated by red star) and the bridge(Br), top 
of boron(TB), hollow(H), top of nitrogen(TN) sites are cal- 
culated as 0.07eV, 0.95eV, l.OOeV and 0.03eV, respectively. 



these two materials, namely carbon atomic chains and 
monolayer h-BN, are combined fundamentally interesting 
properties are attained for promising future applications. 
Interestingly, the properties of the grown structures de- 
pend on the number of carbon atoms in the chains, such 
that they exhibit an even/odd disparity. In addition, 
we showed that CACs grown on h-BN constitute chemi- 
cally active sites for Au, Li and H2, which are normally 
weakly bonded to h-BN. For example, H2 approaching to 
the top of the CAC(2) is dissociated and separated into 
two H atoms, each attached to the chain. We also showed 
that CACs can attach to another h-BN from its free end 
to serve as pillars to increase the spacing between h-BN 
flakes, where molecules like H2 can be stored. We pre- 
sented the electronic energy band structure for various 
lengths of chains to see their variations with the num- 
ber of chain atoms. We finally studied the problem from 
the reverted point of view and presented the growth pat- 
terns, optimized structural and electronic structures of 
BN chains grown on graphene surface. Our results re- 
vealed the unique self-assembly character of carbon and 
BN chains on single layer honeycomb structures and in- 
teresting features attained thereof. 
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FIG. 2: Snapshots of the molecular dynamics simulation 
showing the formation of a short chain comprising four car- 
bon atoms. The snapshots correspond to the initial, 20th, 
40th and 120th steps of the molecular dynamics simulation 
done at 500K. Note that the formation of CAC(4) takes place 
as the CAC (3) leaves its initial bonding position and attaches 
to a single carbon ad-atom at close proximity from its top. 
Similar growth mechanism is also seen during the formation 
of CACs of length n < 8. 



II. METHOD 



In our calculations we use the state-of-the-art first- 
principles plane- wave calculations within the density- 
functional theor}^^ combined with ab-initio, finite 
temperature molecular-dynamics calcu lations using pro- 
jector augmented wave potentials ! 28 * 29 * The exchange cor- 
relation potential is approximated by the generalized gra- 
dient approximation with Van der Waals correction! 30 * 31 * 
Numerical computations have been carried out by us- 
ing VASP software! 3 - 2 ^ A plane- wave basis set with energy 
cutoff of 600eV is used. The Brillouin zone is sampled 
in the k-space within the Monkhorst-Pack schemep^ and 
the convergence of the total energy and magnetic mo- 
ments with respect to the number of k-points are tested. 
The convergence for energy is chosen as 10~ d eV between 
two consecutive steps. In the relaxation of structures 
and band structure calculations, the smearing value for 
all structures is taken as O.OleV. We consider adsorption 
of chains on (4 x 4) supercells and treat the system us- 
ing periodic boundary conditions. The pressure on each 
system was kept smaller than ~ 2kBar per unit cell in 
the calculations. In the ab-initio MD calculations the 
time step was taken as 2.5 fs and atomic velocities were 
renormalized to the temperature set at T = 500K and 
T = lOOOif at every 40 time steps. In the MD stability 
tests, the simulations were run for lOps. 
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III. ADSORPTION OF SINGLE CARBON 
AD- ATOM ON H-BN 



Before going into detailed studies of carbon chains, 
we first investigate the adsorption and migration of sin- 
gle carbon ad-atom, which is the starting point of CAC 
growth on h-BN substrate. We use a (4 x 4) supercell of 
h-BN that consists of 16 boron and 16 nitrogen atoms. 
There is an energy difference of ~ 0.2eV between the spin 
polarized and spin unpolarized energy values in favor of 
the spin polarized state, indicating that the system has 
a magnetic ground state. Therefore all of the calcula- 
tions mentioned hereafter are performed using the spin 
polarized conditions. 

The most favorable binding site of single carbon atom 
was determined by placing the ad-atom initially to var- 
ious adsorption sites at a height of ~ 2A from the 
BN atomic plane and running fully self-consistent geom- 
etry optimization calculations by keeping the ad-atom 
fixed in x— and y— directions and letting the vertical 
z— coordinate of the ad-atom, which is its height from 
the plane, free. Meanwhile, the atoms in the BN super- 
cell are relaxed in all directions except for one corner 
atom of the supercell, which is fixed in all directions to 
prevent h-BN from sliding. 

In Fig. [TJa) the nitrogen and boron atoms of the op- 
timized h-BN structure are separated from each other 
by 1.45 A. The most favorable bonding site of single car- 
bon ad-atom, which turns out to be near the top site of 
nitrogen atom, is marked with a red star. The migra- 
tion(diffusion) path of carbon ad-atom on h-BN with a 
minimum energy barrier is shown by black stars. The 
minimum energy barrier is calculated as 0.68eV^ The 
energy landscape calculated over the whole BN hexagon 
also shows that the energy barrier to the diffusion fur- 
ther increases as the carbon atom moves away from the 
nitrogen atom as shown in Fig. [ljb) . The variation of 
energy calculated along the 2D path shown in Fig. [TJa) 
is presented in Fig. [TJc). As indicated in the figure, the 
most favorable site for the carbon atom is near the top 
nitrogen site, although not exactly on top of nitrogen. 
The energy barrier between the most favorable site and 
the top bridge(BR), top boron(TB), hollow(H) and top 
nitrogen(TN) sites were calculated as 0.07eV, 0.95eV, 
l.OOeV and 0.03eF, respectively. 

In addition to the diffusion path analysis of a single 
carbon ad-atom, we next study the interaction between 
two carbon atoms on h-BN surface. When the distance 
between the two ad-atoms becomes less than a threshold 
distance of ~ 2A, these two carbon atoms attract each 
other and form a CAC (2) perpendicularly attached to h- 
BN. This is indeed the nucleation for longer CACs. The 
most favorable binding site CAC (2), is again near the 
top site of nitrogen. A complete site analysis was also 
performed to confirm this result, by placing a CAC (2) on 
various adsorption sites and comparing the total energy 
values. 



2S\ — r 



- - 



. 2.0 



5? ]5 



0." 

c 

'— 
•- 



1.0 



0.5 



n — \15 




1 2 3 4 5 6 7 & 
Number of carbon atoms in the chain 

FIG. 3: Binding energies (Eb), and the heights(Zi) of odd 
and even numbered CAC(n)'s from the atomic plane of BN 
are shown in green, red and blue lines, respectively. The h 
values exhibit an even/odd family behavior depending on the 
number of carbon atoms in the chain. The sudden peak in the 
binding energy arises from the change of the magnetic state 
of CAC (2) from magnetic to nonmagnetic when it binds to 
hexagonal BN. 



IV. GROWTH OF CARBON ATOMIC CHAINS 
ON H-BN 

A. Chain growth and even/odd disparity 

Growth of the CAC further continues when a third 
carbon ad-atom is introduced at the close proximity of 
CAC (2). However, this time a complete site analysis of 
CAC (3) shows that the most favorable bonding site is the 
top of boron atom, instead of the nitrogen site. The for- 
mation of CAC (3) happens as follow: CAC (2) leaves its 
initial bonding position, moves higher from the BN plane 
and in the mean time it gets closer to the single ad-atom 
until they are bound to each other near the new ener- 
getic site, which is the top of boron site. Similar chain 
growth behavior is also seen during formation of CACs 
at different lengths. This process is further investigated 
with ab-initio molecular dynamics simulations at 500K 
and the snapshots taken from the growth of CAC (4) is 
presented in Fig. |2j We initiate the MD simulation by 
placing a carbon ad- atom and a CAC (3) to their bond- 
ing sites as shown in Fig. |2|i). The simulation was run 
for 2000 time steps and snapshots taken from the initial, 
20 th , 40 th , and 120^ time steps are shown. As the simu- 
lation proceeds to the 2000 t/l step, the chain stays in its 
position shown Fig. [2^iv) , which is an indication of its 
stability at that bonding site. 

We further continue the above analysis with fourth, 
fifth, sixth and seventh carbon ad-atoms. In general, 
once a new carbon atom is implemented at the close prox- 
imity of the existing chain CAC(n — 1), this chain leaves 
its binding site and is bound to the top of new carbon ad- 
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FIG. 4: Side and top views of the most favorable binding configurations of CAC(n) on hexagonal BN are shown in (a) and (b). 
N, B, and C atoms are represented by blue, green and brown balls, respectively. CAC(n)'s having an even number of carbon 
atoms (even n) is bound to BN near the top of nitrogen atom, whereas CAC(n)'s with odd number of carbon atoms (odd n) 
prefer top boron site, with the exception of single carbon ad-atom. The geometries are calculated for a (4 x 4) supercell and 
their stabilities are also tested with MD simulations at T = 500i^ for lOps. In (b), only the carbon atom of the chain that is 
closest to the BN plane is shown. 



TABLE I: Most favorable binding sites, total binding 
energies (Eb), magnetic moments(/x), heights (ft) of CAC(n) 
from the BN plane, and the distances of the lowest car- 
bon atom of the chain from the nitrogen(<ic-Ar) and the 
boron (cLc-b) atom in the BN plane for different n values of 
carbon chains. The bonding sites and magnetic properties of 
CACs on BN exhibit an even/odd disparity. With the ex- 
ception of the single carbon ad-atom, even numbered CACs 
bind to BN near the top of nitrogen (TN) atom and the odd 
numbered CACs bind near the top of boron(TB) atom. Addi- 
tionally, the even and odd numbered chains grown on BN have 
magnetic and nonmagnetic(NM) ground states, respectively, 
with the exception of CAC(l) and CAC(2) cases. All calcu- 
lations were performed on a 4 x 4 supercell which contains 32 
carbon atoms. 



Structure Site E b (eV) /x(/xb) h(A) (fc-iv(A) dc-s(A) 

BN+C ~TN 1.36 2.00 1.50 1.59 1.73 

BN+CAC(2) ~TN 2.19 NM 1.32 1.46 1.83 

BN+CAC(3) ~TB 0.28 NM 1.59 2.16 1.62 

BN+CAC(4) ~TN 0.24 2.00 1.37 1.50 1.63 

BN+CAC(5) ~TB 0.20 NM 1.60 2.28 1.62 

BN+CAC(6) ~TN 0.18 1.97 1.41 1.52 1.58 

BN+CAC(7) ~TB 0.14 NM 1.67 2.41 1.64 

BN+CAC(8) ~TN 0.08 1.96 1.48 1.56 1.59 



atom in their new most favorable binding site. Namely, 
this binding site keeps changing between the near top 
nitrogen site and the near top boron site. With the 
exception of the single carbon ad-atom, we observe an 
even/odd disparity of the binding site depending on the 
number of atoms in the CAC, that is even numbered 
CACs bind to h-BN substrate near the top nitrogen site, 
and the odd numbered CACs bind near the top boron 
site. This situation can be attributed to the charge den- 



sity of CAC(n) at its end, which depends whether n is 
even or odd. 

The binding energies of CACs attached to graphene 
are calculated using the expression = ^[free- 
linear CAC(n)]+£ T [bare h-BN]-£ , T [CAC(n)+h-BN]. In 
this expression, ^[free-linear CAC(n)] is the total en- 
ergy of fully relaxed linear chain of length n, i^x[bare 
h-BN] is the total energy of bare (4x4) supercell of h-BN 
and £?T[CAC(n)+h-BN] is total energy obtained when 
CAC(n) binds to the supercell of h-BN. The complete 
list of these binding energies and the binding geometries 
of these CACs are presented in Table [T[ 

It is seen that, with the exceptional case of CAC(2), 
which is further discussed in the following section, the 
binding energies (E^) of the CACs tend to decrease as 
the length of the chains (or the number of carbon atoms) 
increase. Similarly, due to the decrease in the binding 
energy, the height (ft), namely the distance between h-BN 
plane and the first carbon atom of the CAC, increases. 
However, the ft values also show an even/odd disparity. 
In other words, there is an increasing trend in the ft values 
when the CACs are grouped as even and odd numbered 
chains, but the height of CAC(2n) is always less than 
CAC(2n — 1), as shown in Fig. [3] This situation can 
be explained by the fact thats the C-N bond length is 
shorter than the C-B bond length, and CAC(2n) is a 
more stable structure since it has a more symmetrical 
charge distribution as compared to CAC(2n — 1). These 
decrease ft distances of the even numbered CACs that 
are bound to the nitrogen atom. 

We finally test the stabilities of these h-BN+CAC(n) 
complexes at finite temperatures with molecular dynam- 
ics simulations. The ab-initio MD calculations were car- 
ried out for lOps at T = 500K. The final stable bonding 
configurations of various CACs are shown in Fig. [4] After 
MD simulations, these CACs on h-BN remain stable and 
they don't detach from the BN plane, but they slightly 
pull the nearest B and N atoms of h-BN upwards from the 
plane, as seen from the side views given in Fig.Qa). Also, 
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the chains are not perfectly perpendicular to the BN 
plane, but are slightly tilted at different angles. Fig.[4^b) 
shows the top views of the CACs, where only the carbon 
atoms closest to the BN planes are shown. As noted 
above, the even/odd disparity in the bonding sites of the 
chains with the exception of single carbon atom can be 
seen. It should be noted that ab-initio MD simulations 
carried out for lOps is not sufficiently long to accumulate 
enough statistics. However, calculations at high temper- 
atures are worthwhile to give a hint about the stabil- 
ity at relatively lower temperatures. An unstable chain 
would easily breakdown or detach from substrate. This 
did not happen in this case. Additionally, calculations 
using conjugate gradient method resulted in chain struc- 
tures bound to the substrate with a significant binding 
energies, especially for short CACs. Hence, these tests 
provide sufficient evidence that the chains we studied are 
stable and remain attached to the single layer h-BN. 
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FIG. 5: Variation of C-C bond lengths in the short carbon 
atomic chains CAC(n) grown on BN. While the even num- 
bered chains (blue) have alternating short and long bonds, 
the odd chains (red) have a symmetrical bond length variation 
around their central atom. Note that the number of bonds in 
a chain is n — 1 



B. Electronic and magnetic properties 

Having found the structural properties of CACs on h- 
BN, we next focus on the variations of magnetic and 
electronic properties of h-BN with various lengths of 
CACs attached. First, we perform both spin polarized 
and spin unpolarized energy minimization calculations 
for free (bare) CACs without attaching them to h-BN 
and compare the minimum energies of magnetic and non- 
magnetic cases. In a similar manner with the optimized 
structures of CACs attached to h-BN, it turns out the 
magnetization of free CACs also depends on the num- 
ber n of carbon atoms. Namely that a free chain has 
magnetic ground state if it has even number of carbon 
atoms, and a non-magnetic ground state if it has odd 
number of carbon atoms with the exception of single 
carbon atom in vacuum. Such an even/odd disparity 
in the magnetic moments of finite size CACs was previ- 
ously reported. 19 Hence, when these chains are adsorbed 
to h-BN sheet, the overall magnetic properties of the h- 
BN+CAC(n) structures are maintained with the excep- 
tion of CAC(2). While all other CACs preserve their 
magnetic ground states when bonded to h-BN, CAC(2) 
flips from the magnetic ground state to the non-magnetic 
ground state. This change in the magnetic moment of 
CAC(2) increases its bond strength with h-BN and causes 
a higher binding energy. This is the reason for the sudden 
jump in Fig. [3] 

The disparities observed so far are closely related to 
the bond length variation of CACs with even and odd 
number of atoms First, we note that the final stable 
structures of short CACs grown on h-BN have nonuni- 
form bond length distribution. While the even numbered 
chains have polyyne type of structure with alternating 
short and long bonds, the odd numbered chains have 
symmetrical bond length variations around their center 
points. In other words, the bond length variations for 
the odd chains reach to zero at the center, whereas the 



variation of bond lengths continues in chains having even 
n. This situation is related to the geometrical linear sym- 
metry of the chain as also observed in a recent study for 
bare chains.^ For the odd chains, the geometrical center 
of the chain is the middle carbon atom since there are 
even number of bonds. On the other hand, even chains 
have odd number of bonds and the geometrical center is 
on the middle bond. Since these chains are bonded to 
h-BN substrate only from one side they are free to move 
in the other direction. Hence, the chains need to be sym- 
metrical around the center, and symmetrically equiva- 
lent bonds must have the same length. As a result, the 
two equivalent bonds around the center atom of the odd 
chains have the same length. Variations of bond lengths 
of h-BN+CAC(n) for n=5-8 are presented in Fig. [5] 

This overall symmetry in the geometrical structure of 
the odd numbered chains also leads to a symmetry in 
the distribution of the magnetic moments, making the 
ground states singlet for odd chains. In the singlet state, 
magnetic moments exerted by the electrons cancel each 
other leading to zero net spin and zero magnetic mo- 
ment. Therefore, these structures have a nonmagnetic 
ground state. On the other hand, even chains have an 
uneven distribution of magnetic moments on atoms lead- 
ing them to have magnetic ground states. This magnetic 
preference also alters the binding sites of the chains ac- 
cordingly. One would expect this situation to be reverted 
for chains saturated from both ends because the bound- 
ary conditions for the chains are altered as shown in a 
recent studyP^ 

An individual CAC(n) attached to h-BN can modify 
the electronic band structure of h-BN. If an adsorbed 
CAC(n) is sufficiently far from others, it gives rise to lo- 
calized states in the band gaps and resonance states in 
the band continua of h-BN. In this study we consider a 
single CAC(n) adsorbed to h-BN using a model where 
one CAC(n) is attached to each repeating (4 x 4) super- 
cell with the condition that interactions between adjacent 
CACs are negligible. This model recovers approximately 
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FIG. 6: (a) Electronic energy band structures of CACs grown on h-BN calculated for n— 1, 2, 3 and 4. In the magnetic cases, 
spin up and spin down bands are represented by blue and green lines, respectively. The bands below the red dash-dotted lines 
are fully occupied. The localized impurity states appearing as flat bands in the band gap originate from the p bands of the 
carbon atoms that are at the edges of chains, (b) Isosurfaces of the difference charge densities of chains where yellow and green 
regions designate charge accumulation and charge depletion, respectively. The isosurface values are taken as 0.01 electron/A 3 
for C, C 2 , C 3 and as 0.005 electron/A 3 for C 4 . 



the (4 x 4) folded bands of bare h-BN, except that the en- 
ergy difference between the top of the valance bands and 
the bottom of the conduction bands gradually increases 
from 4.5eV to A.SeV. Additionally, flat bands due to 
CAC(n) occurs in the band gap. These flat bands actu- 
ally corresponds to the localized states due to CAC(n). 
In Fig. [6Fa) we present the electronic energy structure of 
h-BN+CAC(n) calculated for n=l,2,3 and 4 using super- 
cell geometry explained above. The energy positions of 
the filled and empty flat bands in the gap are closely re- 
lated with the energy levels of the corresponding CAC(n), 
which vary with n. For magnetic n=l case, spin- up states 
(i.e. flat bands) originating mainly from C-p xy orbitals 
occur above the top of valence band. Empty spin-down 
states of the same orbital character occur near mid gap. 
In the non-magnetic case of n=2 spin-degenerate C-p xy 
and C-p z states are filled and occur near the top of the 
valence band. While a filled C-p xyz state touches the 
top of valence band both for non-magnetic n=3 case and 
magnetic n=4 case, the positions of empty C-p xy states 
in the band gap strongly depend on n. For all n's the 
resonance states occur in the valence band within the 
energy of 2eV from the top. 

The difference charge density is calculated by subtract- 
ing the charge densities of bare h-BN and bare CAC(n) 



from that CAC(n)+h-BN system by keeping the atomic 
configuration of adsorbed and bare CAC unaltered. In 
Fig. [6^b) the isosurfaces of difference charge density indi- 
cate that important changes from the corresponding bare 
CAC is found where CAC is bound to h-BN. 



C. Applications of CACs on BN 

Apart from its interesting and fundamental aspects, 
the growth of CACs on single layer graphene is a reality 
and they can modify graphene's physical and chemical 
properties. Experimental observations and images ob- 
tained using high resolution TEM 25 indicates that CACs 
can be formed on monolayer graphene surfaces. In this 
study, we show that BN is also a suitable substrate for the 
growth of CAC(n). Pre vious experimental studies 4 and 
density functional resultd^^ have shown that epitaxial 
graphene grown on boron nitride surface is possible by 
means of chemical vapor deposition technique or expo- 
sure of BN layer to carbon atoms or dimers at long time 
intervals. In these studies, it was noted that defected 
graphene structures with undesired chains can also form 
during the growth of graphene ! 20 * 22 * Thus, since the tem- 
perature range of our MD calculations are suitable for 
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FIG. 7: Functionalization of BN sheets through adsorption of carbon chains. For example, a CAC(2), which is strongly bound 
to h-BN, creates chemically active sites for Au, Li and H atoms. H 2 molecule approaching to CAC(2) from sides dissociated to 
form two C-H bonds, whereas O2 remains totally inactive. Ti atom takes the carbon atoms with itself and forms TiH2- 




FIG. 8: CAC(2) and CAC(3) grown between two BN flakes. 
The optimized spacing between the flakes increase from 
3.lA to 4.34A and 5.82A upon the formation of chains. 



atomic layer deposition techniques, it is conjectured that 
similar atomic layer deposition or molecular beam epi- 
taxy methods can be used for chain growth on BN layers 
at specific ambient conditions. 

BN is in general not a chemically active material, how- 
ever with the inclusion of certain ad-atoms its chemical 
activity can be improved. CACs grown on BN can be 
used for functionalization of h-BN layers which will en- 
able the adsorption of certain ad-atoms like hydrogen, 
lithium and gold. This functionalization will especially 
be useful where the contact on an insulating surface like 
BN is desired along with connection of electrodes. In or- 
der to increase the chemical activity of the inactive BN 
layer we use CAC(2) chain. In the previous sections, 
we showed that CAC(2) binds to h-BN with a very high 
binding energy as compared to other CACs. Besides be- 
ing bonded strongly to h-BN, it also creates a chemically 
active sites for the absorption of other ad-atoms. 

Here we test the adsorption of Au, Li, Ti, H 2 and O2 on 
CAC(2) and show that the chemical activity of BN can 
be enhanced through CAC(2) attached to it. Ad-atoms 
were placed in the vicinity of CAC(2) on h-BN and fully 
self consistent calculations were performed. Before intro- 
ducing the ad-atoms, the CAC(2) was placed freely in 
its most favorable position as calculated in the previous 
section. It was seen that Au and Li atoms move towards 
the chain and are bound to the chemically active site of 
CAC(2) with binding energies of 1.77 and 1.53eF, respec- 
tively. These are significantly higher than the binding en- 
ergies of Au and Li on bare h-BN, which are calculated 



as 0.01 and 0.12eF, respectively. This is an important 
result indicating that gold electrodes can easily be con- 
nected to semiconducting h-BN through CACs. A hy- 
drogen molecule approaching to CAC(2) from the sides 
moves upwards and dissociates to form two C-H bonds. 
The binding energy for H atoms is calculated as 3.62eV. 
However, when O2 is introduced to the system instead 
of H2, it stays completely inactive and stay away from 
CAC(2). Finally, we consider Ti ad- atom. Normally, Ti 
is bound to bare h-BN with an energy of 0.80eV\ How- 
ever, when h-BN is functionalized with CAC(2), although 
Ti ad-atom initially is bound to the carbon atoms, it 
doesn't stay there but takes away these CAC(s) from the 
BN layer and forms a TiC2 structure which moves away 
from BN. This final configuration is energetically more 
favorable than Ti binding to h-BN by 5.9eV. The final 
relaxed geometries of all these structures are presented 
in Fig. [7| 



Pillared graphene structures and CACs passivated by 
graphene surfaces from both ends have also been subjects 
of recent interest of both theoretical and experimental 
studies.^ It was recently shown that CACs passivated 
by graphene flakes from both ends produce highly stable 
chain structures.^ Oxidized graphene pillars were also 
shown to be useful materials for storage applications™ 



Motivated by previous studies, we also calculated the 
bonding geometries of CACs between BN layers and 
showed that CACs can also be grown between two BN 
flakes as shown in Fig. [8j We have demonstrated this sit- 
uation by calculating optimized bonding configurations 
and the spacing between two BN flakes when CAC(2) 
and CAC(3) are grown between them. CAC(2) is bonded 
to nitrogen atom from one side and to boron atom on the 
other side. On the other hand, CAC(3) is bonded to both 
of the flakes from the top of boron atom. Once CAC(2) 
and CAC(3) are grown, the spacing between the flakes 
increases from 3.1 A to 4. 34 A and 5. 82 A, respectively. 
The formation energies, calculated by subtracting the en- 
ergy of two BN planes and the energy of CAC from the 
energy of the final structure, are found as 2.0 and 0.32eV 
for CAC (2) and CAC (3), respectively. 



TABLE II: Most favorable binding sites, binding energies (i£b) 
and the heights(Zi) of B, N ad- atoms and BN or NB atomic 
chains from the graphene plane. Here, for example, BN(3) 
indicates the BN chain consisting of three atoms (namely N- 
B-N) grown on graphene with N being attached to graphene. 



Structure Binding site Eb(eV) h(A) 



B-ad-atom 


Bridge 


1.24 


1.45 


N- ad- atom 


Bridge 


1.05 


1.25 


NB(2) 


Bridge 


2.58 


1.46 


BN(2) 


Top 


2.42 


1.32 


NB(3) 


Bridge 


1.18 


1.52 


BN(3) 


Top 


0.90 


1.31 


NB(4) 


Bridge 


0.82 


1.58 


BN(4) 


Top 


0.46 


1.39 


NB(5) 


Bridge 


0.73 


1.63 


BN(5) 


Top 


0.21 


1.48 


NB(6) 


Bridge 


0.51 


1.67 


BN(6) 


Top 


0.15 


1.57 






BN ATOMIC CHAINS GROWN ON 
GRAPHENE 



FIG. 9: Growth dynamics of short NB chains on graphene. 
Once a new ad-atom is placed at a random point in the neigh- 
borhood of a NB chain having n atoms, the structure rear- 
ranges itself to form a chain comprising n + 1 atoms. The 
snapshots are taken from the quantum MD simulations per- 
formed at 500K for lOps. Blue, green and brown balls are, 
respectively, nitrogen, boron and carbon atoms. 



Having shown the stable carbon chain structures grown 
on BN, we next investigate the growth of short BN chains 
on graphene to see if a similar self-assembly mechanism is 
also present for this reverted situation. Here we adopt the 
representation where BN(n)+graphene indicates a BN 
chain consisting of n atoms that is attached to graphene 
through N atom. For example, NB(3) + graphene (or 
NB(3) only) stands for B-N-B chain attached to graphene 
through B atom. Following the similar procedures de- 
scribed in Sections III and IV, we calculate most favor- 
able binding sites, electronic properties, binding energies. 
We first perform MD simulations to see the self assem- 
bly mechanisms of BN chains. To this end we start by 
placing single B and N ad-atoms on graphene supercell 
and run MD simulations at 500K for lOps. Eventually, 
the two ad-atoms form a BN molecule perpendicularly 
attached to graphene. When a third ad-atom N is placed 
on the graphene at the close proximity of existing BN 
molecule, MD simulations show that a chain comprising 
three atoms, B-N-B forms. This procedure was repeated 
for different lengths and it was observed that B and N 
atoms are indeed self- assembled on graphene honeycomb 
structure to form short BN atomic chains. Snapshots 
taken from these MD simulations are presented in Fig. [9] 

The most favorable binding sites of short BN chains on 
graphene and their corresponding binding energies were 
calculated with self consistent conjugate gradient calcu- 
lations. It was found that, the binding site of a BN chain 
on the graphene depends on the type of its closest atom 
to the graphene. A NB(n) chain is bound to the graphene 
layer at the bridge site, and a BN(n) is bound at the top 



site as shown in Fig. 10 However, both single B and sin- 
gle N ad- atoms are bound at the bridge site. Contrary 
to the CAC case, no clear even/odd disparity was ob- 
tained in h-BN+graphene complexes owing to the ionic 
bonding between B-N atoms. Here there is another ef- 
fect related with the number of electrons in the chain. 
Chains having odd number of atoms have odd number of 
electrons, while even numbered chains have even number 
of electrons. The binding energies of these chains shown 
in in Table [H] suggest that BN chains are bound with 
higher energies to the graphene as compared to the bind- 
ing energies of CACs on BN, especially for longer chains. 
The binding energy of a boron ending chain (i.e. NB) 
is always higher than a nitrogen ending chain (i.e. BN) 
site. Here we note that the supercell size and adsorbate- 
adsorbate 

The even/odd disparity observed in the magnetic prop- 
erties of CAC(n)+h-BN are also not seen for BN chains 
grown on graphene. Single B ad-atom as well as sin- 
gle N ad-atom adsorbed on graphene have magnetic 
ground states. However, while BN(2) and NB(2) chains 
on graphene are nonmagnetic, BN(3) and NB(3) on 
graphene have magnetic ground states. This order is 
disturbed by the chain consisting of four atoms; while 
BN(4)+graphene is nonmagnetic, NB(4)+graphene has 
magnetic ground state. 

The electronic structures calculated for a single ad- 
sorbate, namely single boron and nitrogen ad- atoms, 
and single NB(2) and BN(2) molecules, perpendicularly 
attached to (4x4) supercell of graphene are shown in 
Fig. 11 Owing to their odd number of electrons, sin- 



9 




FIG. 10: The most favorable binding configurations of short NB and BN chains grown on graphene. A NB(n) chain is attached 
to graphene at the bridge site, since its nearest atom to graphene is B, and BN(n) is bound at the top of carbon atom, since 
its nearest atom to graphene is N. Energies are calculated using a (4x4) supercell and their stabilities are tested with MD 
simulations at T = 500i^ for 10 ps. Note that the most favorable binding site is the bridge site for both single N and single B 
ad-atoms. Blue, green and brown balls are, respectively, nitrogen, boron and carbon atoms. 



gle N and B ad- atoms shift the Fermi level up and down 
relative to the Dirac point of graphene. This is an im- 
portant feature to be used in electronic applications of 
single layer graphene. For NB(2)+graphene the bands 
crossing at the Fermi level attribute semimetallic charac- 
ter. On the other hand, BN(2)+graphene is a semicon- 
ductor with a narrow indirect band gap. This situation 
shows the dramatic effect of the type of chain atom which 
is attached to graphene. In addition, the supercell size, 
adsorbates-adsorbate coupling, as well as the symmetry 
of the adsorbate+supercell system are crucial f or the re- 
sulting electronic structure for small supercells! 41 * 42 * For 
example, the semiconductor presented in Fig. [IT] (d) be- 
comes metallic when BN+graphene system is treated in 
(6x6) supercell. In the case of very large supercell having 
single adsorbate attached the band crossing is recovered 
and the states of adsorbate appear as localized (reso- 
nance) impurity states. 



VI. CONCLUSION 

In conclusion, using first-principles calculations within 
the density functional theory, we revealed the growth 
mechanism of carbon and BN atomic chains on single 
layer honeycomb structures, namely hexagonal BN(or h- 
BN) and graphene, respectively. We found that with the 
inclusion of each new carbon atom, the existing atomic 
chain consisting of n carbon atom leaves its previous po- 



sition to join with the new carbon atom to from a chain 
of n + 1 carbon atom. Similar, but more complex growth 
mechanisms were also found for BN chains grown on 
graphene layers depending on the type of chain atom(B 
or N) which is attached to graphene. These growth pro- 
cesses were simulated by ab-initio molecular dynamics 
calculations at various temperatures and the resulting 
chain structures were shown to be stable even though 
they are free to bend and slightly tilt. Nevertheless, these 
simulations are performed at high temperatures(lOOOK) 
and the chains did not detach from the substrate plane for 
lOps which is adequately long for ab-initio molecular dy- 
namics simulations. Therefore these simulations presents 
reliable stability results, especially at room temperature. 
The growth of atomic chains on the single layer honey- 
comb structures heralds a self-assembly process, which 
may have fundamental and technological implications. 

The grown chains by themselves, exhibit interesting 
physical and chemical properties depending on the num- 
ber of atoms forming the chain and the type of chain atom 
attaching to the substrate. In particular, the physical 
properties of even numbered and odd numbered carbon 
chains can behave differently leading to interesting even- 
odd disparity. We also showed that atomic chains grown 
on single layer substrates attribute useful functionalities 
to bare h-BN and graphene. These properties are depen- 
dent on the number of carbon atoms in the chains. Apart 
from creating localized electronic states in the band gap 
and local magnetic moments it is demonstrated that car- 
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(a) B-Graphene (b) N-Graphene 




(c)NB-Graphene (d)BN-Graphene 
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bon chains can also be used for increasing the interlayer 
spacing between BN flakes, where specific molecules can 
be stored. Atomic chains grown on graphene and h-BN 
create chemically active sites on BN surface for atoms 
such as Au, Li and H 2 , which provides the connection 
of these materials. Additionally, specific atoms forming 
strong bonds with atomic chains modify the electronic 
structure in the band gap and hence change the conduc- 
tivity for possible sensor applications. Similarly, modifi- 
cations of the electronic properties of graphene depend- 
ing on the grown B or N ad-atoms and BN or NB chains 
are also worth emphasizing. Moving the Fermi level of 
graphene up or down by adsorbing N or B atoms, re- 
spectively can be a worthwhile feature in the electronic 
applications related with graphene. The self-assembly in 
terms of atomic chains grown on semimetallic graphene 
and semiconducting h-BN and interesting functionalities 
achieved thereof can bring up new perspectives for fur- 
ther research on single layer honeycomb structures. 
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